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ASaTRACT 

We  investigate  the  excitation  and  propagation  of  the  three  dimensional 
electromagnetic  field  over  an  infinite  corrugated  plane  which  is  approximated 
by  an  anisotropic  impedance  boundary  condition.  Emphasis  is  placed  upon 
effects  of  surface  anisotropy  which  are  not  evident  in  two  dimensional  treat¬ 
ments,  In  particular  we  consider  the  excitation  by  a  magnetic  point  dipole 
in  detail.  It  turns  out  that  the  fields  are  determined  by  a  scalar  wave 
function  which  satisfies  a  mixed  boundary  condition  involving  a  linear  com¬ 
bination  of  the  wave  function,  its  normal  derivative  and  its  second  order 
tangential  derivative.  The  exact  formal  solution  is  first  derived,  and  then 
the  radiated  far  field  and  the  surface  wave  far  field  are  evaluated 
separately.  Both  the  phase  and  the  an5>litude  of  the  excited  surface  wave 
are  dependent  upon  the  direction  of  observation.  Numerical  results  are 
given.  The  physical  significance  of  this  solution  is  discussed.  A  com¬ 
parison  is  made  between  this  problem  and  the  theory  of  ship  waves. 


1.  Introduction 

In  this  paper  we  study  the  excitation  and  propagation  of  the  three 
dimensional  electromagnetic  field  over  an  infinite  corrugated  plane  which 
may  be  approximated  by  an  anisotropic  impedance  boundary  condition.  This 
investigation  is  motivated  by  its  potential  application  to  scanning  surface 
wave  antennas^.  The  fields  produced  by  a  magnetic  point  dipole  will  be 
considered  in  detail.  The  boundary  condition  is  such  that  surface  waves 
are  generated  and  it  may  be  described  mathematically  as: 
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one  involving  only  a  single  Fourier  transform  by  introducing  an  auxiliary 
function  which  is  similar  to  those  recently  employed  for  several  other 

li  R  6  7 

problems  in  diffraction  theory  ’  ’  ’  . 

In  section  2  we  obtain  the  formal  solution  for  the  fields  of  a 
magnetic  point  dipole  oriented  in  the  direction  of  corrugations  and  located 
above  a  corrugated  plane  which  is  approximated  by  an  anisotropic  impedance 
boundary  condition.  The  dipole  is  first  decomposed  into  phased  line  soxurces, 
and  the  fields  of  the  line  sources  are  then  synthesized  into  the  dipole 
fields.  In  section  3  the  far-zone  radiation  field  is  evaluated  by  an  algebraic 
method  involving  asyitptotlc  differentiation  and  the  excited  surface  wave  far  field 
is  determined  by  applying  the  method  of  stationary  phase.  All  the  surface 
wave  characteristics  which  include  the  aii5)litude,  the  phase,  and  the  ex¬ 
ponential  decay  factor  are  functions  of  the  direction  of  observation  and  the 
surface  reactance.  In  section  4  we  summarize  the  peculiar  properties  of  the 
three  dimensional  surface  wave,  and  present  the  surface  wave  power  patterns 
and  the  surface  wave  phase  patterns.  A  comparison  is  also  made  between  this 
problem  and  the  theory  of  ship  waves. 

2.  The  Formal  Solution 

In  this  section  we  proceed  to  study  the  electromagnetic  field  which 
arises  from  a  magnetic  point  dipole  lying  above  the  surface  of  a  plane 
characterized  by  an  euiisotroplc  impedance  boundary  condition.  Consider  a 
magnetic  point  dipole  located  at  the  point  X  =  0,  y  =  h,  z  =  0  (See  Fig.  l) 
and  oriented  in  the  direction  of  the  z-axis.  We  are  Interested  in  the  free 
space  region  y  >  0.  The  Impedance  boiuidary  condition  prescribed  on  the  plane 
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surface  y  =  0  is  given  by 

(2.1)  E  =  0  E  =  ZH=-iXH 

Z  X  z  z 

where  the  impedance  Z  is  supposed  to  be  purely  inductive  and  will  support 
surface  waves.  We  nov,'  wish  to  solve  the  time  reduced  Maxwell's  equations 
subject  to  the  prescribed  singularity  and  boundary  conditions. 

If  we  assume  the  time  dependence  to  be  of  the  form  e”  ,  the  mono¬ 
chromatic  Maxwell's  equations  may  be  reduced  to: 

^2  2)  curl  H  =  —  iwe  E 

^  curl  E  =  icJMH-5(x)6(y-h)6(z)z 


where  €  and  M  are  the  permittivity  and  magnetic  permeability  of  free  space. 
Because  of  the  boundary  condition,  the  z-conrponent  of  the  electric  vector 
is  absent  in  the  fields  and  hence  the  field  is  completely  determined  by  the 
value  of  a  scalar  wave  function  f.  We  have 
(2.3)  E  =  curl  f  z  and  H  =  curl  otu:l  f  z 


i.e. , 
(2.4) 


E 

z 


0,  E 


M 

by 


5x 


and 


(2.5) 


H 


H 


1 

icon 

bx.bz 

1 

icon 

bydz 

1 

icon 

Ldz2 

The  function  f  satisfies  the  equation 
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(2.6) 


(v^  +k^)f  =  6(x)  6(y  -  h)  B(z) 


where  is  the  three-dimensional  rectangular  Laplacian,  k  is  the 
propagation  constant  of  free  space,  and  8  is  the  ULrac  Delta  function. 

Ihe  appropriate  boundary  condition  maybe  deduced  from  (2.1)  and  (2.3): 

2 

(2.7)  +  —  (^^^+1:1)  =  0  when  y  =  0 

ay  uH  dz2 

X 

where  —  is  a  positive  number.  Bie  mathematical  problem  has  been  reduced 

to  that  of  solving  the  inhomogeneous  wave  equation  (2.6)  subject  to  the 

mixed  second  order  boundary  condition  given  by  (2.7).  Furthermore,  the  far 

field  should  satisfy  the  radiation  condition. 

Applying  Fourier  transform  to  eqs.  (2.6)  and  (2.7),  we  obtain 
2  2 

(2.8)  (\+  +  K^)  F  =  8(x)  6(y-h) 

8y 


(2.9) 


^  +  2L  k2  F  =  0 


where  the  functions  f(x,y, z)  and  F(x,y,s)  are  related  by  the  following  Fourier 
transform  pair: 


(2,10) 


f(x,y,z)  =  f  F(x,y,s)  e^ 

-  00 

00 

(2.11)  F(x,y,s)  =  J'  f(x,y,z) 


s  z 


ds 


-i  s  z  , 
e  d  z 


p  2  2 

and  IC  =  k  -  s  .  Mathematically  the  transformed  equations  are  a  two- 
dimensional  inhomogeneous  wave  equation  and  a  mixed  first-order  boundary  con¬ 
dition.  Physically  the  point  dipole  has  been  decomposed  into  phased  line 
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sources.  We  are  going  to  first  solve  the  problem  of  a  phased  line  source 
and  then  synthesize  the  line  source  fields  into  the  dipole  field. 

The  field  of  a  phased  line  sovirce  may  be  analyzed  either  by  the 
Fourier  transform  method  or  by  using  another  auxiliary  function  which  is  a 
linear  combination  of  the  unknown  wave  function  and  its  derivative.  We 
choose  to  employ  the  latter  method  because  it  gives  another  exan^ile  of  this 
relatively  new  technlg.ue.  Let  us  introduce  the  auxiliary  function, 

N 

(2.12)  G  =  +  X)  F 

oy 

where  X  =  —  K^.  Then  G  satisfies  the  equation 
to  n 

2  2 

(2.13)  (-^  +  -^+K^)G  =  -  (—  -  X)  B(x)  5(y-h) 

dx  by  5h 

and  the  simple  homogeneous  boundary  condition 
(2.  Ilf)  G  =  0  when  y  =  0 

2  2 

We  first  treat  the  case  of  positive  X,  i.e.,  s  <  k  .  If  the  function 
F  is  denoted  by  F^  when  X  is  positive,  then  F^  can  be  found  from  G  by  Integrating 
(2.12)  in  the  following  manner. 

y 

(2.15)  F^(x,y,x)  =  ®  ^  ^  y'  G(x,T1,s)  dy  +  e“^^  T(x,s) 

—  00 

where  T(x,s)  is  an  arbitrary  function  of  x  and  S,  which  except  for  certain 
discontinuities  in  the  derivative  must  be  a  homogeneous  solution  of  (2.I3). 

The  discontinuities  must  be  adjusted  so  that  F^(x,y, s)  meets  the  requirements 
of  the  problem.  This  adjustment  follows  after  eq.  (2.I8)  below. 
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By  applying  the  operator  X)  to  the  lonique  Green's  function 
satisfying  eq.  (2.8)  and  vanishing  at  y  =  0,  we  obtain  the  solution  for 
the  auxiliary  function: 


(2.16)  G  =  f  (^  -  X) 

4  oh 


0^’  [kn/xS*  (y-h)^] 


,H<1> 


^'/x^  +  (yfh)^ 


We  substitute  (2.l6)  into  (2.I5)  and  introduce  the  excited  surface 
wave  terms  to  play  the  role  of  e" 

y 


T(x,s) 
(l)t 


(2.17)  F^(x,y,s)  =  -I  ^(Tl.h)2]  +  |j  ^  [yx2+(Tl+h  )2] 


fd  T1 


i 


tf 

I 


>11 


(1)  , _ , 

Ho  LVx2^<Tl-h)2J 


-  H 


(1) 


LVx^  +  (Tl+h)^. 


+  C  e'^y  Wk2  +x2  |x| 

Integrating  by  parts,  we  have 

(2.18)  F^(x,y,»)  .  -  J  -  i  nf’-)  [xy^-riT^l^] 

e’'"  H<^>  [ltV;^7nPr^2]  dT| 

-09 


+  C  g-Xy  +  d  /  U  I 
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Now  the  surface  wave  amplitude  C  may  he  determined  hy  requiring  the 
^1 

continuity  of  -r —  across  the  line  x  =  0,  i.e., 
ox 


(2.19) 


■  Sli  1 


Lim 

x-»0 


iF^(x>0,  y).  y)j=  0 


Substituting  (2.l8)  into  the  Jun®  condition  (2,19)  yields 

y 


(2.20) 


dTl 


+  12  C  v'k2  +  3^2  ^\/k^  +  |; 


The  limit  of  the  integral  in  the  above  equation  contributes  onl>  in  the 
neighborhood  of  T|  =  -  h.  Hence 


(2.21) 


X  K  e 


•Xh 


-h  +e  ( 1) 


C  = 


2y  k2  + 


Lim  r 
X-+0  J 


H 


1  Tv^^  +  cn+h)^]: 


+  (T1  +  h)^ 


:dTl 


-h-  e 

Using  the  asymptotic  formula  of  Hankel  function  for  small  argument,  we  find 
that 


-h  +  e 


(2.22)  C  = 


,  rr  _-X  h  n  -21  xdT| 

=  ■  ■  Lim  /  ( —  ■=)  -==== 

2Jk^  +  x2  x-»0'^  511^x2  +  (Tl+h)2  ^x^  +  CH  +  h 


-h-  € 


Carrying  out  the  Integration  gives 


iX  e 


-Xh 


(2.23) 


C 
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2  2 

Now  we  consider  the  case  of  negative  X,  i.e.,  s  >  k  .  If  the 

function  F  is  denoted  by  F^  when  X  is  negative,  then  Fg  may  be  found  from 

G  by  Integrating  (2.12)  as  follows: 

00 

(2.24)  V  F2(x,y,s)  =  -  f  G  (x,Tl,s)  dT] 

y 

Substltxrtlng  (2.l6)  into  (2.24)  and  integrating  by  parts,  we  have 

(2.25)  F2(x,y,s)  =  ^[Ky^2  +  (y.^)2]  -  |  ^  [k^^2  ^  ^ 

e^^  [k  ^^2  +  (T)  +  h)^]  IT) 

y 

The  limits  of  integration  in  (2.24)  have  been  so  chosen  that  they 
insure  the  convergence  of  the  integral.  Because  _ i.  is  continuous  across  the 

8  X 

the  line  x  =  0,  no  siirface  waves  are  excited  in  this  case.  Since  the  arguments 
of  Hankel  Sluictions  in  (2.28)  are  imaginary,  Fg  will  not  give  rise  to  a 
radiation  field. 

The  substitution  of  (2.l8)  and  (2.25)  into  (2.10)  completes  the  formal 
solution  of  this  problem.  It  1s  difficult  to  evaluate  the  integrals  exactly^ 
however,  we  shall  obtain  the  far  field  expressions  in  the  next  section. 

3.  Determination  of  the  Ihr  Fields 

In  this  section  we  determine  the  radiated  far  field  and  the  surface 
wave  far  field.  We  first  treat  the  radiated  far  field  by  an  algebraic  method 
involving  asymptotic  differentiation.  This  method  has  proved  useful  in  two 


f 
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dimensional  prdblems^'^'^’'^.  It  gives  the  radiated  far  field  directly, 
without  intermediary  use  of  the  previously  obtained  complete  solution. 
Let  us  now  introduce  the  following  function 


(3.1)  g 

Ihen  g  satisfies  the  equation 


^f  X  .  2  S2f  , 

5-y  *  ^  *  ^72) 


(3.2)  (V^  +k^)  g  =  [j  ^  ^2^]  &(y-h)  6(2) 

and  vanishes  when  y  =  0.  The  solution  for  g  may  be  written  down  immediately: 

(3.3)  g  =  ^ 


'  A  ^ 

i  (k=  .  4 

r  ikR’  ikR"n 

IF".-®  1 

-  dh 

W|i  d  z‘ 

J  L  r'  r”  J 

where  R'  =  'Jj?  +  (y  -  h)^+  «  R  -  h  sin  0  sin  ^  , 

R' '  =\/x^  +  (y  +  h)^  +  K  R  +  h  sin  0  sin  (|)  , 

R  =  \/x2  +  y^  +  z^  , 


1  2  1  y* 

0  =  tan"  ■  ,  and  ^  s--  tan"  - 

v/x2+  y^  X 


We  may  expect  the  feir-zone  radiation  field  of  f  to  be  of  the  form 


(3.4) 


f  = 


M0^)  ikR 


Substituting  (3.3)  and  (3,4)  into  (3. l), introducing  polar  coordinates,  and 
keeping  terms  of  order  we  obtain 


{[• 


[-Ik  sin  0  sin  <)  -  Iji  (kV  cos^  0)Je"^^  ®  ^ 

-I  ikR 

(3.5)  t  [-tt  .in  8  .l„  f  +  JL  (fcS.kS  .In  8  sin  ^  — 

^  k  sin  0  Bln  ^  (ii^  -  oos^  0 )  J  f 
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(3.6)  f 


gikR  p  -ikh  sin  0  sin  ij)  .  ^  ^  -  n  sin  6  g^jj  q 

— —  e  '  +  '  '"■ '  ■  —  ■  ■  ■  e 

^  Jt  R  I-  i  sin  (j)  +  n  sin  0 


where  n 


^  X 

“v/I7e 


may  he  interpreted  as  normalized  surface  reactance. 


Now  we  t\zm  our  attention  to  the  surface  wave  far  field.  Substituting 
the  last  term  of  (2.l8)  into  (2.10)  and  using  (2.23),  we  obtain  the  surface 
wave  field; 


(3.7)  fg 


r  ^  .-Xy  +  i  [v/k^  +  |x|  +  s  z]  ^  ^ 

\  \Iy^  +x2 

-  k 


It  should  be  remembered  that  X 


—  and  =  k^ 


2 

s  . 


In  order 


to  evaluate  the  integral  by  the  method  of  stationary  phase  it  is  convenient 
to  make  the  following  substitutions 


(3.8)  z  =  p  cos  |x|  =  P  sin  f  =  Q! 


(3.7)  becomes 
(3.9)  fg  = 


1 

J  kif,  n,  a)  e^  ^  ^ 
-1 


P(t,  n,  a) 

d  P 


where 

(3.10) 


A(^,  n,  a) 


2  p 

ink  e"“  ^(^  "  ^  ^  (^ 

2Jt  \jl  -  +  n2(l-a2)2 


and 

(3.11)  P  (.f,n,a)  =  a  cos  f 


sin  f 
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when  k  p  is  sufficiently  large 


(3.12)  fg  =  A(ag) 


r - 2« - 1 

|J!  P(V',n,a3  )  I  -1 


^  i  [k  p  P(V',n,a  )"  i  ] 


It  can  te  easily  shown  that  - x  P(^,n,ag)  Is  always  negative.  The 

Sa 

stationary  point  a  mist  satisfy  the  following  equation: 
s 


A  p  etc  1 1  +  2a 

(3.13)  —  =  cos  ^  -  sin  ^ 


r  2  2  1 

1  +  2n  (l-gg) 


da 


[  1  -  a^  +  n^  (i-q{  )2  ji 


=  0 


We  may  easily  express  V'  in  terms  of  n  and  a 


(3.14)  tan  f 


1  -  ag^  +  n^  (l-Og  f 

s  . . . ii.i.gfii 


a  [1  +  2  n^  (l-a  J  ] 


The  right  hand  side  of  equation  (3.14)  is  a  monotonically  decreasing 
function  in  the  interval  0  <  ct  <  1  for  all  n.  The  truth  of  this  statement 
is  suggested  hy  Fig.  2  and  can  be  proved  by  differentiating  (3.14)  with  respect 
to  a  .  Therefore  a  can  be  found  numerically  for  the  directions  of  observation 

S  6 

0  <  ^  <  90°.  Those  values  corresponding  to  the  interval  90°  <  ^  <  180°  and 
-1  <  a  <0  are  obtained  by  the  symmetry  requirement.  Now  the  surfhee  wave  feu: 
field  may  be  written  as 


g-nk  (l-Q^)  (y  +  h) 
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where  a  Is  a  function  of  ij/  as  indicated  by  (3,14).  The  electric  and  magnetic 
s 

field  components  of  the  far-zone  surface  wave  may  be  obtained  by  substituting 

(3.15)  into  (2.4)  and  (2.5).  However,  in  order  to  understand  the  power  flow 
of  the  surface  wave  far  field,  it  is  desirable  to  use  the  cylindrical  co¬ 
ordinates  (P,i>}y).  We  first  write  down  the  components  of  f  : 

8 

,  ,  ,  A  A  A 

(3.16)  f  z  =  f  COS  ^  p  -  f  sin  ^  ^ 

S  S  6 

Substituting  (3.16)  into  (2.3)  yields 

(3.17)  4*  =  P  [-n  k(l-a^)  sln^f  J  +  ^  f-  n  k  (l-CK^)  cos  ^  f  J 

+  y  r  -  1  k  P  sin  ^  f  -  i  k  cos  ^  f  1 


(3.18)  H  '  p  [k^  sin  f  +  k^  (||  )^  cos  f  f^  -  n^L^d-a^)^  cos  f  f 

+  f  [n^^  (l-a^)^  sin  f  -  k^P^  sin  f  -  k^P^^cos  f  f  J 
+  y  r  -  i  n  k^(l-oP)P  cos  ^  f  +  i  n  k^(l-Q!^)  ^  sin  ^  f  1  , 

S  S  S  ^’jh  s 


1 

where  all  terms  of  higher  order  than  p”^  have  been  neglected.  The  average 
Poynting  vector  is 


(3.19)  iaeCEiH  )  .iae(ly 


"y  -  ®y  >V  )  0 


+  i  Re  (Ey  Hp*  -  Ep  Hy*)  ^  i  Re  (Ep  E^*)  y 
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Lengthy  but  straightforward  substitutions  into  the  above  equation  show  that 


the  f  and  y  components  vanish  and  give  the  expression  for  the  p  con^xjnent 


(3.20).  (P  )„  =  i 

'  av  p 


(l-oj)  [  [l  +  2  n^  (i-a^)  ]^  [i-a^  +  n^d-a^)^]  • 


where  f  is  given  in  (3.15). 
s 


Integrating  (3.20)  with  respect  to  y  gives  the  pattern  function  for  the 
total  surface  wave  power. 


exp  [-2nk  (l  —  Q^)hJ 


This  is  a  function  of  f  through  equation  (3.14). 

4.  Discussions 

Both  the  phase  eind  anplltude  of  the  excited  surface  wave  found  in  the 
last  section  are  dependent  upon  the  direction  of  observation.  The  siirface  wave 
amplitude  shows  the  typical  exponential  decay  away  from  the  reactive  surface. 
Notice,  however,  that  the  exponential  decay  factor  also  depends  upon  the  di¬ 
rection  of  observation.  Because  the  direction  of  observation  is  approximately 
constant  in  any  fixed  small  region  far  away  from  the  source,  the  three  di¬ 
mensional  surface  wave  may  be  regarded  locally  there  as  an  oblique  plane  surface 
wave.  It  is  recalled  that  none  of  the  anisotropic  characteristics  of  the  dipole- 
excited  surface  wave  are  present  in  the  two  dimensional  treatments  where  an 
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infinite  line  source  is  assumed  ’  or  where  the  propagation  of  an  oblique 
plane  surface  wave  is  considered^.  In  the  case  of  a  finite  line  source  we 
may  expect  to  observe  the  two  dimensional  behavior  at  points  close  to  the 
line  source;  however,  the  three  dimensional  behavior  should  appear  when  the 
distance  from  the  source  is  large  compared  with  the  length  of  the  source. 

Ihere  must  be  a  transition  between  the  two  dimensional  and  three  dimensional 
characterl stlc  s . 

The  surface  wave  phase  pattern  and  the  surface  wave  power  pattern 
have  been  computed  numerically  for  different  values  of  normalized  surface 
reactance  n  and  h  =  0.  The  phase  curves  using  rectangular  coordinates  are 
plotted  in  Figure  3  for  0  <  ^  <  90°.  Those  corresponding  to  90°<  f  <  180° 
may  be  obtained  by  the  symmetry  requirement.  The  phase  curve  for  n  =  0.1 
is  very  close  to  unity.  The  surface  wave  power  patterns  using  polar  coordinates 
are  plotted  in  Figures  4  and  5  Tor  a  magnetic  point  dipole  of  unit  strength. 

The  siuface  wave  power  density  vanishes  at  ^  =  0  and  180°  for  all  values  of  n. 

It  is  interesting  to  notice  that  the  surface  wave  power  patterns  for  n  =  0. 1  and 
0.5  are  broadly  directional  with  the  maximum  at  f  =  90°,  while  those  for  n  =  5 
and  10  are  sharply  directional  with  the  maximum  at  an  angle  close  to  ^  =  0° 
and  180°.  This  observation  may  lead  to  the  conclusion  that  the  surface  wave 
power  pattern  of  a  finite  line  source  resembles  that  of  a  continuous  array  of 
isotropic  point  sources  for  small  values  of  n.  However,  a  similar  statement 
is  not.  valid  in  the  case  of  large  values  of  n. 

If  frequency  is  very  low,  the  limiting  form  of  the  wave  equation  and 
boundary  condition  in  our  problem  will  be  reduced  to  the  aame  as  those  in 
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Kelvin's  ship  wave  prohlem.  But  we  are  interested  in  the  fields  above 
the  reactive  plane,  while  the  ship  wave  solution  deals  with  the  space  below 
the  water  surface.  Although  the  wave  amplitude  and  phase  in  both  problems 
depend  upon  the  angle  of  observation,  we  find  one  wave  in  all  directions 
but  vanishing  a.t  f  =  0^  and  l80°,  while  the  ship  wave  consists  of  two 
orthogonal  waves  and  is  confined  within  a  sector  behind  the  ship. 
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Figure  4.  Surface  Wave  Power  Angular  Density  x  (Sttcu/z.) 


Figure  5.  Surface  Wave  Power  Angular  Density  xlSTTW/i). 
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